clear all; close all
global nx0 ny0 nz0 dif Lwin amp Lper tau pert dt Ls 
tic
amp=1/5;
restart=0;             
nts_start=0; 
tau=5e-7;

ntp=500; nts=800;dt=ntp*tau; 
Ls=1;Lper=5;
pert=1e-3; 
alx=1;nx0=50;
%aly:y方向长度;ny0:y方向格点数
aly=5;ny0=250;
alz=20;nz0=32;
Lwin=alx/8;eta=1;
saldir =sprintf(['./data/']);
saldir1 =sprintf(['./figure/']);
parameters;
scmod;
% visualize;
visual1;
toc
% showt;